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FORWARD 



This report serves as an introduction to the subject of human 
describing function measurement. Particular emphasis has been placed 
upon development of the spectral analysis relations utilized in the 
describing function measurement techniques. The work was performed by 
Dr. Hess as part of a research study sponsored by the Air Force Flight 
Dynamics Laboratory, Air Force Systems Command . 

This technical memorandum has been reviewed and is approved. 
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ABSTRACT 



A review of the spectral analysis techniques used in the 
measurement of human describing functions is presented. The 
describing function relations for single loop, compensatory tracking 
tasks are developed. The effect of sinusoidal inputs and finite run 
lengths are discussed and a brief discussion on mechanization 
techniques is included. 
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I, INTRODUCTION 



This memo is a brief introduction to the subject of human 
describing function measurement with particular emphasis placed upon 
mathematical formulation. The presenta’. ion is intended to serve as a 
guide to the uninitiated. 

The problem to be treated here is one of modeling the human 
operator via the describing function technique in a single axis, 
compensatory, tracking task. For discussion of such tasks, see Ref. 1, 
pp. 7-9. Figure 1 is a block diagram of the compensatory system. 
Experiments have shown that in tasks such as this, the human operator 
is nonlinear and time varying in behavior. However, a good deal of 
success has been achieved in treating him in a quasi-linear fashion 1 . 
This quasi-linear model assumes that, for the most part, the human 
behaves in a linear, time invariant manner. This means that the major 
portion of his response can be attributed to a linear, time invariant 
operation on his visual stimulus. 

Figure 2 is a more detailed block diagram in which the operator 

1 2 R 

has been modeled via a linear describing function and a remnant * y . 
The term describing function is preferred to transfer function to 
emphasize the fact that the model is approximating a nonlinear element 
and is valid only for particular inputs. Nonetheless, the term transfer 
function does appear in the literature; e.g.. Ref. 4. 

Briefly then, the describing function accounts for the operator’s 
linear behavior, the remnant, his nonlinear time varying behavior. 

Being more specific, in Figure 2: 

p(t) represents that portion of the total operator output 

c ( t ) which can be obtained by a linear operation on the 



1 



System 

Forcing Operator Operator System 




FIGURE 1 

The Human Operator in a Compensatory Tracking Control Task 
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Nonlinear Operator 








FIGURE 2 

Equivalent Block Diagram of Human Operator 
in Compensatory Tracking Control Task 
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visual stimulus e(t). 

n(t), the remnant, represents that portion of the total 
output c(t), which is not linearly correlated with 
the input i(t). 

It should be evident that the utility of this quasi-linear 
technique depends upon the extent to which the operator is, indeed, 
linear. If the remnant is relatively large, the describing function 
technique is of little value in itself. 

The modeling problem to be treated here is an empirical one; i.e., 
to determine (j(u) and §^(co) based upon the physical measurements 
of finite tracking runs in laboratory experiments. The conditions 
under which these measurements are to be taken are as follows: 

1) The conditions implicit in Figure 2 are in effect; 
i.e., the task is single axis, compensatory, and the 
system dynamics are linear and time invariant. 

2) The input i(t) is truly random or at least random 
appearing. 

3) The operator is well trained. This simply means that 
his adaptation and learning periods have passed. 

Of these three conditions, perhaps the second deserves a brief 
comment. The reason that the input form must be established is that 
the operator's dynamics are a definite function of the input; e.g., 
the operator's behavior vhen tracking a single sinusoid has been found 
to be considerably different than when tracking a random or random 
appearing input. Indeel, the very nature of the task can change when 
the input is predictable . ^ The rationale behind choosing a random or 
random appearing input is that it more truly represents the environment 
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i„ kmn operator models are to be applied; e.g., pilot pitch 

attitude tracking in the presence of atmospheric turbulence. , 
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II. SPECTRAL ANALYSIS 



The techniques used in describing function measurement depend 
to a great extent upon spectral or harmonic analysis . The brief 
review which follows is intended to summarize these analytic tools. 

For a more thorough treatment the reader is referred to Ref. 6 5 
Chapters 2 and 13. 

A. Periodic Signals 
1. Fourier Series 

A periodic signal x(t)* with fundamental frequency 
which satisfies the Dirichlet conditions (Ref. 7, p, 248) can be 
represented by a Fourier series 

<x> 

x(t) = ^ X(n)e J " n<B l fc . (l) 

n=-® 

T/2 

X(n) = | J dt (2) 

-T/2 



T 



2rr 

®L 



2. Autocorrelation and power spectal density 
The autocorrelation function for the periodic signal x(t) is 
defined as 



T/2 

V (t) = |J x ( t ) x ( t+T ) dt (3) 

-T/2 
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It can be shown that 



(>*) 



»„(") - |x(n)| 2 



(5) 



Where § (n) is referred to as the power spectral, density of the 



xx 



signal x(t). It can be shown that 



T/2 

W») - W d * 

- T/2 



( 6 ) 



and from eqn. (4) 



=L W n ^ T (7) 

n=-oo 

3. Crosscorrelation and cross power spectral density 

The crosscorrelation function for two periodic signals 
x(t) and y(t) with identical fundamental frequencies and which 
satisfy the Dirichlet conditions is defined as 

T/2 

^(t) = | J x(t)y(t+-r)dt (8) 

-T/2 

It can be shown that 



^(t) = ^ X(n)Y(n)e' 1 " nu) l T (9) 

n=-oo 
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X(n) denotes the complex conjugate of X(n). Kow, 






8 (n) = X(n)Y(n) (10) 

Where ^(n) is referred to as the cross power spectral density of 
the signals x{t) and y(t). One can show 

T/2 

§ (n) = J f cp (r)e~^ n \ r dr (11) 

xy T J xy 

-T/2 

and from eqn. (9) 

CO 

v ( t) = Z V n)ednv (12) 

Two periodic signals are said to be linearly uncorrelated when 

cp (t) = 0 for all t. 

xy 

B. Transient Signals 

1. Fourier Integral 

A signal x(t) is said to be transient if 
lim x(t) = 0 

If such a signal 

1) satisfies the Dirichlet conditions in any finite Interval 

2) satisfies the inequality 

CO 

j |x(t)|dt < OO 

_co 

Then the signal can be expressed as a Fourier integral (Ref. 7, p. 279) 
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(13) 



x(t) = J X(ouj)e' 3 *' t ' t dto 



where 



X(j(u) = j x(t)e^^ at 



(14) 



2. Autocorrelation and energy spectral density 

The autocorrelation function for the transient signal 
x(i) is defined as 



cp^Ct) = J x(t)x(t+T)dt 



(15) 



It can be shown that 



<Pxx ( t) = h J ! x(jU)) i 2e0 J>T du> 



(16) 



Now, 



*^(<0) = |x(d«>)i : 



(17) 



Where §^.(u)) is referred to as the energy spectral density of the 



signal x(t). It can be shown that 

co 

§ xx (lu) = J c fkx (T)e_DU)T dT 



(18) 



and from eqn. (l6) 



c Pxx ( t) = h J § xx ((B)eJUJT dtu 



(19) 
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3. Crosscorrelation and cross energy density spectra 

The crosscorrelation function for too transient signals 
x(t) and y(t) each of which satisfies the Dirichlet conditions in 
every finite interval and which satisfy 

CO 

J |x(t)j dt < « 

-CO 

CO 

J |y(t)J dt < ® 



is defined as : 



CO 

^(t) = J x(t)y(t+T)dt 
-00 



( 20 ) 



It can be shown that 

CO 

^(t) = ~ J X(ju))Y(jto)e jU5T do (21) 

-co 



Now 



= X( jto)Y(ju)) 



( 22 ) 



Where § (w) is referred to as the cross energy spectral density 

xy 

of the signals x(t) and y(t). 

It can be shown that 



V> * I 



(23) 



and from eqn. (21) 
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CO 



v ( t) = h! v (cu)sJU)T du) {2h) 

— CO 

Two transient signals are said to be .early uncorr elated when 

cp (t) = 0 for all t. 



C • Random Signals 

1. Fourier Integral 

Consider a random signal x(t) as a sample function from 

a stationary* ergodic random process. Since* in general, 

00 

I" jx(t)|dt 



is not finite* one cannot write a Fourier integral representation 

for x(t). 

2. Autocorrelation and Power Spectral Density 

The autocorrelation function for the random signal x(t) 



is defined 



T 

^(0 = lira It J x ( t ) x ( t+T ) dt 



( 25 ) 



Now can be represented by a Fourier integral since it 

satisfies the too conditions of section II. B. 1. Hence* it can be 
shown 

oo 

*XX (t) = h J § X3C (U5)eJCOT dU> (26) 

— CO 
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where $ (cd) is referred to as the power spectral density of the 

XX 

signal x(t). Here 



§ 

xx 




-jun* 



dr 



( 27 ) 



3 . Crosscorrelation and Cross Power Spectral Density 

The crosscorrelation function of two signals x(t) and 
y(t) which are sample functions from two different random processes* 
each of which are stationary and ergodic and jointly ergodic* is 

defined 



T 

<p (t) = lim §T J x ( t )y( t+T ) dt 

** T-* 00 m 

-T 



(28) 



Now it can be shown that 

00 

c p (t) = ~~ f § (u^e^ 7 dtu (29) 

•xy 2rr J xy 

.00 



where §^(u>) re ^ erre( i t o as the cross power spectral density of 
the two random signals x(t) and y(t), Here 



§ (u>) — f cp (r)e dr 

xy J ^xy v ' 



( 30 ) 



Two random signals are said to be linearly uncorrelated when 



cp (t) = 0 for all r. 
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0. Ergodicity 



Briefly, ergodicity is the property of a random process which 



ensures the equivalence of statistical and time averages. For example, 
if x(t) is a sample function from an ergodic random process 



i.e., the expected value of the signal x(t) at time (a statistical 
average) equals the time average on the right hand side of the 
equation. 



measures (they’re now equivalent) taken from a single experiment 
are indeed representative of the random process under study. 

E. Relations for Linear Systems 

Ref. 6, p. 333 shows that, given a linear, time -invariant 
system, with weighting function (impulse response) h(t), input x(t), 
and output y(t), where x(t) and y(t) are the random functions alluded 
to in section II. C. 2, 3 



T 




Ergodicity also ensures that the statistical or spectral 



x(t) v>- h(t) £>- y(t) 




( 31 ) 
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Alt 



H(joj) 



5 («>) 

xy 

$ ((I)) 

xx 



(32) 



where 



CO 

H(j(u) = J h(t)e~^ U ^ > dt 

— CO 



CO 



J h(t)e~ Jl£ * dt = H(s)j 
i I s=g”<2> 



( 33 ) 



and H(s) is the system transfer function. 
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m. DESCRIBING FUNCTION RELATIONS 

A. Problem Statement 

Consider Figure 2 again. Let the input i{t) he a sample 
function from an ergodic random pro. ;s. On the basis of measurements 
of finite time histories of the signals i(t), e(t), c(t) and m(t), 
one wishes to find 

1) The operators linear describing function Y (ju>). 

P 

2) The power spectral density of the remnant signal n(t), 

§ (u>). 

nn 

B. Finite Run Length 

Obviously, for the purposes of practical data collection, 
only finite time histories of the system signals of Figure 2 are 
available. But referring to section II. C., one sees that integration 
over infinite time is implied in the various spectral measures of the 
random signals. The question of just how long the histories must be 
in order to ensure accurate spectral measurements naturally arises . 
Here it will simply be stated that T, the run length, must be long 
enough to ensure accurate determination of the autocorrelation 
functions . 

Consider the system input, i(t), in Figure 2. Define 



i T (t) = i(t) - T < t < T 



i T (t) = 0 



other t 



Now let 



'P ii ( T ) T = §jT j i T (t)i T (t+T)dt 

-T 



( 34 ) 
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Not that 

s 



lim cp..(t) t = <p (r) 

*p— ICO 



Also define 



’ f dT 



-T 



Note that 



(35) 



(36) 



lim § (<u) = § (w) 

T-co 11 T 



(37) 



In light of the definition of i T (t), eqns . (34) and ( 36 ) can he written 



*ii (T) T = It J V t}1 T (t+T)dt (38) 



§ ii M T = J Cp ii (T) T e "° t ° T dT 



(39) 



Finally , since i,p(t) looks very much like a transient signal, one can 
assume that the conditions of section II. B. 1 are applicable and that 
the Fourier integral (or transform) of i^(t) can be written; i.e., 

CO 

1 ( 3 ®) = J i T (t)e~ jwt dt (40) 
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Furthermore, from eqns. ( 38 ) and (39) 

CD CO 

$ ii (t0) T = J 2T J i T ( t ) i T ( t+T ) e " J&T dtdT (*H) 



-CD — 00 



Now defining a new variable u as 



u = t + T 



one can write 



$ ii^T = It J / i T ( t ) i T (u)e" Jtl,i ' U “ t ^ dtdu (42) 

-CD _ CD 

- rl ! *1 / *»] 



= jbj£ i ^ U) ^ j J = f^liU^)! 2 (^3) 



from eqns. ( 37 ) and (43) 



§ii(u)) = lim[ |^|l(ju))| 2 j 



m 



•Likewise, one can show 



$ ic (u)) = lim[ l( ju))C ( ju>) j 



(^5) 



J J' — *CD 
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The limiting processes of eqns* (44) and (45) should be 
interpreted as allowing T to be large enough to ensure accurate 
spectral measurements but finite so that eqn. (4o) is still valid, 
C. Spectral Relations 

Referring again to the block diagram of Figure 2 



or 



E(ju>) = l(jo)) - M(jo?) 

= I(joj) - [ N(jo>) + E(ju>)Y (>) ] Y c (>) 



I(ju)) - N(jiu)Y (jcu) 

E(jco) = 



1 + Y (ou>)Y (jto) 

c p 



Now 



l(jcu)E(ju)) = 



I(jiu) - N(jto)Y c (j(u) "j 

1 + Y (ju>)Y (ju>) 

P 



In like manner , dropping the (jco) notation,, one can show 

tPyI+Yn! 

p c J 



IC = 



1+YY 
c p 



EC = 



EE = 



f I - Y H 1 


Y I + N I 


L c Jl 


_ c J 


r 1 + y y 


T 1 + Y Y 1 


L c p . 


JL c p J 


r 1 - 7 » 1 


I-Yll 


L c Jl 


c J 



[ 1 + Vp 1 1 * Y = Y p ] 



(46) 



(47) 



(48) 



(49) 



18 



